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The equations for calculating diffraction profiles for bent crystals are revisited 
for both meridional and sagittal bending. Two approximated methods for comput¬ 
ing diffraction profiles are treated: multilamellar and Penning-Polder. A common 
treatment of crystal anisotropy is included in these models. The formulation pre¬ 
sented is implemented into the XOP package, completing and updating the crystal 
module that simulates diffraction profiles for perfect, mosaic and now distorted 
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crystals by elastic bending. 


1. Introduction 

The motivation of this work is the availability of easy-to-use 
computer codes for simulating diffraction or reflection profiles 
of bent crystals. The target audience is the users and researchers 
of synchrotron radiation facilities, but also other fields of X-ray 
research, like plasma physics. There is a vast literature on the 
theory and applicability of the Dynamical Theory of Diffraction 
to the calculation of diffraction profiles (see e.g, (Authier, 2001) 
for an updated review). However, many scientists and engi¬ 
neers need to know the performances of bent crystals under 
X-rays without becoming specialist in the field of Dynami¬ 
cal Diffraction. Several available codes may be used, many of 
them available via collaborations with their authors: REFLECT 
(Etelaniemi et al, 1989), REFLEX (Caciuffo et al, 1990), 
PEPO (Schulze & Chapman, 1995), DIXI (Holzer et al, 1998), 
and others publically available (Stepanov, 2004). One popular 
code for the synchrotron community is XOP (Sanchez del Rio 
& Dejus, 2011), a graphical environment for computer codes 
for i) modeling of X-ray sources (e.g., synchrotron radiation 
sources, such as undulators and wigglers), ii) calculating char¬ 
acteristics of optical elements (mirrors, filters, crystals, multi¬ 
layers, etc.), and iii) multipurpose data visualizations and analy¬ 
ses. XOP is used extensively to simulate crystal diffraction pro¬ 
files for perfect, bent, and mosaic crystals. 

The calculation of diffraction profiles for flat (undistorted) 
crystals is usually performed using the basic equations of the 
Dynamical Theory, in different formulations. Although a full 
quantum treatment of X-ray diffraction by a crystal exists 
(Ashkin & Kuriyama, 1966; Kuriyama, 1967), a classical 
approach, as discussed in (Hartwig, 2001) is adopted for most 
practical cases. In the classical approach, the propagation of the 
X-ray field inside and outside the crystal is described by the 
Maxwell equations assuming that all quantities (electric sus¬ 
ceptibility, electric field, etc.) are defined in a continuous way 
for any point in the space and time, and the interaction of the 
X-ray field with the crystal electrons is described by quan¬ 
tum mechanics. The dynamical theories for X-ray diffraction 
in perfect crystals have been extended to include lattice dis¬ 
tortions related to crystalline defects and macroscopic bending. 


For practical purposes, several levels of approximations may be 
defined (Hartwig, 2001): 

(a) dynamical diffraction theory for the perfect crystal, 

(b) local applications of the dynamical theory for perfect crys¬ 
tals to distorted ones, 

(c) geometrical optics or eikonal theory, and 

(d) wave optics or Takagi theory. 

Item (a) is related to the classical dynamical diffraction of 
X-rays in perfect crystals. The term dynamical is used when 
the rescattering and absorption of the X-rays in the crystal 
volume is considered. Perfect crystals are ideal undistorted 
monocrystals over a large distances (as compared with the 
unit cells) thus assuming perfect alignment of the atoms in 
the crystalline structure. There is no curvature of the atomic 
planes (like the originated by elastic bending or thermal dis¬ 
tortion) and no alteration of the atomic order (no inclusions, 
dislocations, defects, stress, cracks, etc). Since the pioneering 
work of (Darwin, 1914), several several formulations are avail¬ 
able (Ewald, 1917), (von Laue, 1931), (Zachariasen, 1967), 
etc. Some comprehensive books describe them (James, 1994), 
(Pinsker, 1978) and (Authier, 2001), which includes a complete 
historical review. 

The (b) approximation works for crystals with small dis¬ 
tortions with respect to the perfect crystals. This means that 
the displacement vector u{r) varies very slowly. The crystal 
deformation transforms a point with position r into another at 
r', with uif) = r' — r K, u{r')). The crystal reflectivity for 
diffracted beam can therefore be computed quantitatively by 
shifting the angular position of the incident direction by a value 
equal to the effective misorientation, and applying the diffracted 
intensity of the perfect crystal. This method was introduced 
by (Bonse, 1958) and (Authier, 1966). Also, the multilamel¬ 
lar method (Caciuffo et al, 1990; Erola et al, 1990) discussed 
later belongs to this category. It can be used for estimating the 
diffraction profile in both Bragg and Laue geometries in many 
practical cases of crystal curvature. 




The diffraction theory of the (c) approximation mimics the 
geometrical optics for visible light. An inductive derivation was 
made by (Penning & Polder, 1961a), a deductive theory was 
presented in (Kato, 1963; Kato, 1964a; Kato, 1964^), and a 
derivation from the Takagi-Taupin Equations is in (Indenbom & 
Chukhovskii, 1971). Level (d) rely on the Takagi-Taupin equa¬ 
tions (Takagi, 1969; Taupin, 1964a). 

For calculations of perfect undistorted crystals XOP imple¬ 
ments the equations of the dynamical theory of diffraction from 
(Zachariasen, 1967), summarized in Section 2.1. On the other 
hand, imperfect mosaic crystals can be simulated in XOP also 
using Zachariasen theory which is valid for mosaicity values 
much larger than the Darwin width. 

In this work we extended the XCRYSTAL application in 
XOP to cover perfect crystals with small deformations origi¬ 
nated by elastic bending. Elastic bending produces a crystal dis¬ 
tortion that depends on its anisotropy. Every crystalline material 
is anisotropic, meaning the the elastic constants are not scalar, 
but depend on the direction (they are deduced from the com¬ 
pliance or stiffness tensors). Anisotropy is important when the 
crystal is bent, and is not relevant for the perfect undistorted 
crystals. We distinguish curvature in two directions with respect 
to the direction of the incident beam: meridional curvature in a 
plane that contains the incident direction, and sagittal, a plane 
perpendicular to it. The meridional plane coincides with the 
diffraction plane. A cylindrically bent crystal is curved in one 
single plane. A single moment is sufficient to bent the crystal 
in one direction. Two moments along perpendicular planes will 
bend the crystal in two directions. Perfect cylinders are diffi¬ 
cult to obtain by elastic bending, because when applying one- 
moment to a crystal block or plate there is an spurious curva¬ 
ture in the plane perpendicular to the main bending plane (anti- 
clastic curvature). Spherical or toroidal crystals are curved in 
both planes, usually made by applying two moments. The elas¬ 
tic constants and tensors are related to the principal crystallo¬ 
graphic directions, which are coincident with the crystal block 
directions only if the crystals is symmetric, i.e., the crystalline 
planes are parallel to the crystal faces. In the most general case, 
the crystal planes are not parallel to the crystallographic direc¬ 
tions and the crystal is called asymmetric. For the diffraction 
effects, a crystal curved with a non constant radius of curva¬ 
ture (parabolic, ellipsoidal, conic, etc.) can be approximated as 
a crystal with two averaged curvatures over the meridional and 
sagittal planes. 

For simulating the diffraction profiles, two approximated 
theories are used: the multilamellar method and the Penning- 
Polder theory. 

The multilamellar method described in Section 2.2 can be 
used to simulate diffraction profiles of curved crystals in both 
Bragg and Laue geometries. This models belongs to level of 
approximation (b). We develop the formulation of the multil¬ 
amellar theory working for both sagittal and meridional curva¬ 
tures in both Bragg and Laue geometries. This requires a correct 
treatment of the crystal anisotropy. The unified formulation pre¬ 
sented here extends the cases treated in literature for isotropic 
Bragg crystals like (Caciuffo et ai, 1990; Erola et al, 1990) and 


for anisotropic Laue crystals (Shi, 2011) to the general case of 
two-moment bending Laue or Bragg anisotropic crystals. 

The Penning-Polder model (Penning &. Polder, 1961a) sum¬ 
marized in Section 2.3 belongs to the level of approximations 
(c), and only applies to Laue geometry. This model has been 
successfully applied for calculating the diffraction profiles of 
high energy monochromators used in beamlines at many syn¬ 
chrotrons (APS, NSLS, ESRF, Spring-8, Petra, Diamond, etc.). 
Again, we present here a formulation that can be applied for 
crystals bent in two directions (meridional and sagittaly, thus 
unifying previous uses of the Penning-Polder theory for merid¬ 
ionally bent crystals, isotropic (Sanchez del Rio et al., 1997) or 
anisotropic (Schulze &. Chapman, 1995), or for sagittal bend¬ 
ing (Shi, 2011). Approximation level (d), thus solving Takagi- 
Taupin equations will be addressed in a future work. 

From the computer point of view, a unification of and mod¬ 
ernization of the XOP crystal module has been done in order 
to upgrade, clean and improve its structure. A single Fortran 
95 module calculates now crystal diffraction using the different 
calculation algorithms described here, using a full 3D vectorial 
calculus for beam direction and crystal orientation. Thus, this 
paper is a good companion of the software package and will be 
used as reference manual. It is designed for being integrated in 
other X-ray codes, like for the ray tracing package SHADOW 
(Sanchez del Rio et al., 2011). 

2. Algorithms for computing diffraction by bent 
crystals 

Using a crystal reflection defined by the Miller indices hkl, the 
reciprocal vector of the lattice is H = (l/dhki)n^, with dhu 
the interplanar distance, and a unitary vector normal to the 
Bragg planes [hkl]. 

The Laue equation 

k^=kl+H, ( 1 ) 

which is satisfied only for the diffraction condition, gives the 
wavevector of the diffracted wave kf for a particular position 
of the incident wavevector kg, that is, its angle with the reflect¬ 
ing [hkl] planes is the Bragg angle 9b. This gives the Bragg law 
A = 2dhki sin 9b, with A the photon wavelength. In this paper 
|A:| = 1/A as in the text of (Zachariasen, 1967). 

The change in the direction of any monochromatic beam (not 
necessarily satisfying the diffraction condition or Laue equa¬ 
tion) diffracted by a crystal (Laue or Bragg) can be calculated 
using i) elastic scattering in the diffraction process: 

|k0| = |k«| = i (2) 

with Ic'^ = (1/A)y®'^, and V a unitary vector; and ii) the 
boundary conditions at the crystal surface: 

= (3) 

where 11 refers to the component parallel to the crystal surface. 

A crystal cut is defined by n, a unity vector normal to the 
crystal surface pointing outside the crystal bulk. Usually it is 
expressed as a function of a in Bragg geometry, and x in Laue 
geometry, but if both values are well-defined they can be used 



indistinctly in both geometries (see Appendix A). The projec¬ 
tions of the beam directions onto this vector are: 70 = n • V® 
and 7 /r = n • . The asymmetry factor is h = in ■ kP)Hn ■ 
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2.1. The perfect crystal 


For the perfect (undistorted) crystal we follow 
the (Zachariasen, 1967) formulation, because of its accuracy, 
compactness and easy numerical implementation. This formu¬ 
lation expresses the crystal reflectivity as a function of angular 
parameter: 
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which measures the separation of the incident field from the 
Bragg condition either in angular terms (“rotating crystal”, for 
a fixed photon wavelength): 


az « 2(05 - 6) sin(26lB) = -2A6» sin(26lB), (5) 


or in terms of photon wavelength (or energy): (“Laue method”) 
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( 6 ) 


However, it is recommended to compute az using the exact 
expression in Eq. (4) which is also valid in extreme cases, like 
in normal incidence. 

The X-ray reflectivity of a single perfect parallel-sided crys¬ 
tal in Bragg (or reflection) geometry is: 
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For Laue (or transmission) geometry we have: 
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The transmitivity (forward diffracted beam) are 
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( 10 ) 


/*’ is the intensity of the incident wave along direction V® with 
wavevector is the intensity of the external diffracted wave 
along the direction with wavevector ,c 1^2 are phase terms 
dependent on the crystal thickness T, and both X 17 and C 17 
terms depend on az and on the crystal electrical susceptibility 
in a rather non-trivial way: 


q = b'VH'Vfi, P is the polarization factor (P = 1 for a- 
polarization, P = | cos 20^1 for 7 r-polarization), and '^h is the 
Fourier component of the electrical susceptibility 'Fq related to 
the structure factor F)/ as: 
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where ro the classical electron radius, Vc the volume of the unit 
cell, e the charge of the electron and c the speed of light. 

The C 17 phases in Eqs. 7-10 are expressed as: 
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and the other quantities are defined as: 


(13) 


(14) 


Note that a factor |h|“* appears for calculating the diffracted 
beams (Eqs. (7) and ( 8 )) to guarantee the conservation of the 
total power when the linear width of the incident beam is large 
compared with the depth of penetration in the crystal, as dis¬ 
cussed in (Zachariasen, 1967) (pag. 122). This factor is not 
present for the transmitted beams (Eqs. (9) and (10)). Also note 
that the signs of ^12 in Eq. 13 are changed with respect to 
(Zachariasen, 1967) because in our definitions the surface nor¬ 
mal points outside the crystal. 

az in Eq. (4) is the magnitude that measures the separa¬ 
tion of the incident beam k® from the Bragg position. One can 
define a dimensionless parameter 7 that measures the “normal¬ 
ized” angular separation or “deviation parameter” (Zachariasen, 
1967),(Authier,2001). 

vwmi vwmi 

The Bragg uncorrected angle verifies the Bragg law A = 
2dhki sin Ob- The angle at p = 0 corresponds to the Bragg angle 
corrected by refraction 6bc'- 
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Bragg angle and corrected Bragg angle are equal only for Laue 
symmetric case (b = 1). The Darwin width is the angular inter¬ 
val {A6)£) that corresponds to a 2r] width centered at 7 = 0 : 
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The Darwin width corresponds to the total reflection zone 
( 11 ) for a non-absorbing thick Bragg crystal, and to the EWHM 


z = 



(Full Width at Half Maximum) for the non-absorbing Laue 
crystal. It is often used as an indicator of the width of the 
diffraction profile. The Darwin width is exactly the FWHM 
for Laue non-absorbing crystals, but for the Bragg the case 
of thick-nonabsorbing crystals (known as Ewald solution) it 
is FWHM = 2V3/3(A6»)£. « 1.155(A6»)d. (or FWHM = 
3-\/2/4(A0)o « 1.061 (At?)/) for the small-absorbing infinitely- 
thick crystal, or Darwin solution) (Zachariasen, 1967). For a 
general crystal the FWHM can be calculated numerically from 
the simulated diffraction profile. 

Another important parameter is the extinction depth. The 
extinction in crystals is associated to the crystal thickness 
needed to diffract most of the beam, and it is associated to the 
primary extinction coefficient ^ext or attenuation of the incident 
beam due to the diffraction. The extinction depth A is the depth 
at which the intensity of the incident wave is attenuated by a 
factor He: 
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Its value is usually given at the center of the diffraction pro¬ 
file (?7 = 0). In some cases, the extinction is given for amplitude 
instead of intensity, and the value is twice the one defined in 
Eq. (18): A.ampi = 2A. Moreover, one can also define the extinc¬ 
tion length along the incident beam path (instead of depth along 
the crystal normal), thus Aungth = A/| 7 o|. Usually extinction is 
associated to crystals in Bragg geometry, but for weakly absorb¬ 
ing crystals (as for Laue crystals) it is more appropriated to use 
the Pendellosung depth Ape„d {Apend = 27rA). 


The parameter c is a reduced curvature, the “deformation gra¬ 
dient”. It is a constant for uniform bending, and its expression 
comes directly from Eq. (19) (see also (Albertini et al., 1977; 
Taupin, 1964??)): 


c = 


drj 
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( 21 ) 


and depends on the bending geometry and elasticity parameters 
of the crystal. Appendix B gives an introduction to the elastic 
anisotropy in crystals. An general expression of c for a doubly 
curved Laue or Bragg crystal is obtained in Appendix C: 
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where M are the bending moments, Sij are components of the 
compliance tensor (see Appendix B), I is the inertia moment 
of the crystal, and A, coefficients depending on the in and out 
beam directions: 


Al = (y")V3° - (y°)X 

A 2 = y3"y3°(y3" - y3°) (23) 

A3 = y3''y3°(V2'' - Vi) 


2.2. The Multilamellar (ML) method 

The main idea behind this method is to decompose the crys¬ 
tal (in the direction of beam penetration) in several layers of 
a suitable thickness. Each layer behaves as a perfect crystal, 
thus the diffracted and transmitted beams are calculated using 
the dynamical theory for plane crystals. The different layers 
are misaligned one with respect to the others in order to fol¬ 
low the cylindrical surface of the crystal plate. This model was 
first introduced by (White, 1950) and further developed by, 
among others, (Egert & Dachs, 1970) and (Boeuf et al, 1978). 
It has been used for optimization of monochromators for inelas¬ 
tic scattering and coronary angiography applications (P. Suortti, 
P. Pattison and W. Weyrich, 1986; Erola et al., 1990), and 
for simulating crystal analyzers for fusion plasma diagnostics 
(Caciuffo et al, 1990). 

In the ML method, the bent crystal of thickness T is decom¬ 
posed into a series of perfect crystal lamellae of constant thick¬ 
ness AT. The value of rj is then a function of the depth t (t = 
axis) in the crystal from the entrance surface, or 


r]{t) = r]{0) + cA, (19) 


where p(0) is the p value (Eq. (17)) at f = 0 or entrance sur¬ 
face, A is the crystal thickness in units of extinction depth for 
the amplitude given by (Zachariasen, 1967) (Eq. 3.140) 
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( 20 ) 


Based on the adopted model (see Eig. 1), the Bragg planes in 
each lamella are tilted relative to the ones in its neighbor lamella 
by an angle Ap = 7r/2 for the Laue case and Ap — 2 for the 
Bragg case. Therefore, the reduced thickness for the lamella is 
AA = Ap/jcj (Eq. 22), thus giving AA = 7r/(2|c|) for Laue and 
AA = 2/|c| for Bragg case. The thickness of a lamella is (Eq. 20) 
AT = 2A AA and the number of lamellae is = A/AA = T/AT. 

The total reflectivity of the given set of layers can be com¬ 
puted by writing the energy balance for the n layers, which leads 
for the Bragg case to: 

N ( 7-1 
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and for the Laue case: 

TZ = Y,1 rf Y[ t'r^^ \ , (25) 
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where Sh = AT/pi, is the X-ray path of the diffracted beam 
inside a single lamella, p, is the absorption coefficient of the 
crystal material, r, and f, are the reflectivity and transmission 
for the i-th layer, respectively, that are computed using Eqs. (7) 
to (10), depending on the geometry (Bragg or Laue). 






Figure 1 

Crystal division in the Multilamellar model for a Bragg and Laue crystal with 
the respective transmitted and diffracted beams for Bragg (top) and Laue (bot¬ 
tom) geometries. 


In the ML model it is assumed that the beam trajectory inside 
the crystal is a straight line. The crystal is perfect inside a 
given lamella, thus the crystal curvature cannot be large, oth¬ 
erwise it would originate local strains in the crystalline planes. 
The model is valid for crystals sufficiently thick to guarantee 
the existence of several lamellae. The model fails if the crystal 
thickness is of the order of or smaller than the lamella thickness. 
This may occur in Laue cases, where the crystal should be thin 
enough to guarantee a high transmission. The oscillations that 
may be found in Laue profiles calculated with the ML model 
are due to unphysical interferences between the crystal lamel¬ 
lae. A detailed description of the assumptions of this method is 
in (Caciuffo et ai, 1987) and (C.T. Chantler, 1990) and a good 
validation with experiments is in (Erola et al, 1990) 

2.3. The Penning-Polder (PP) method 

The PP theory is another geometrical approach widely used 
because of its simplicity in generating the beam trajectory inside 
the crystals. In the original model (Penning & Polder, 19611?), 
only the two “normal wavefields” are accounted for. It is only 
valid when the disorientation of the crystal lattice planes over 
the Pendellosung length (or extinction length) is much smaller 
than the intrinsic reflection width of the perfect crystal. How¬ 
ever, by including the two “created wavefields” from the inter¬ 
branch scattering, the PP theory can be extended to the case 
of strongly distorted crystals (Balibar et al, 1983; Schulze & 
Chapman, 1995). 

In the ray-optical theory of (Penning & Polder, 19611?) the 


X-ray beam in a distorted crystal is assumed to be a pseudo¬ 
plane Bloch wave (’’wavefleld ray”) propagating parallel to the 
local Poynting vector. The crystal is supposed to be composed 
of parts of flat and undistorted crystals where the dynamical 
theory for perfect crystals can be applied. The wavefleld is pre¬ 
served passing from one part of the crystal to the next. Diffrac¬ 
tion phenomena (the interference between two wavefields) are 
neglected. Thus, the Pendellosung fringes are not simulated 
with this model. 

For crystals under constant strain gradient j3 (see Appendix 
D) through the crystal thickness T, the ratio ^ of the amplitudes 
of the diffracted and the transmitted waves can be obtained from 
Eq. 33 in (Penning & Polder, 19611?) 

277 = e. - bl^e = e- - bl^i + 2/37’, (26) 

where the subscript e stands for the “exit” surface, i denotes the 
“incident” surface. 

The above equation generates two solutions for and another 
two for 

Cc.i = ?? + \/if -b 
ie,2 =21- -b 

e-,1 =V-PT + y'ipT -r^y-b (27) 

6,2 = 2]- /3T - v (/37’ -r]f -b 

these two solutions correspond to the splitting of the wavefields 
at the incident (6;/) 7 = 62) and exit surfaces Most of 
the intensity goes along one mode for plane waves (Balibar 
et ai, 1983). For each mode, the relationships between the 
intensity of the incident beam (/q), the diffracted beam (/«) 
and the transmitted beam {It) are (Eq. 35 in (Penning & 
Polder, 19611?)); 
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The total reflectivity TZ is then obtained by adding the inten¬ 
sities of the two ^ solutions (one is usually very small), and 
removing the intensity of the created wavefields from the inter¬ 
branch: scattering (Balibar et ai, 1983), or 


TZ = TZi + 72.2 




1 — exp 



, (29) 


where /3c = Tr/{2Ape„d) is the critical strain gradient introduced 
by (Authier & Balibar, 1970), with Apend the Pendellosung 
period as defined in Sec. 2.1. In case of overbending (|/3| >> 
/3c), the TZj beams are attenuated and the intensity flows along 
the transmitted beam direction; 
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3. Examples of application 



























3.1. Cylindrically bent crystals in Bragg geometry 


In a first example we checked the output of our code against 
the experimental case described in Fig. 4 of (Erola et ai, 1990), 
where it is shown the reflectivity for symmetrical Si400 reflec¬ 
tion using Mo Kai radiation (E = 17479 eV) for three bend¬ 
ing radii; 1.1, 2.7 and 5.7 m. The bending is in the meridional 
plane, and the crystal is considered isotropic (Poisson’s ratio 
V 2 i = —S 2 \ls 22 = 0.28). The results of our simulations are in 
Fig. 2. The results including crystal anisotropy are almost iden¬ 
tical to the isotropic case. 
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Figure 2 

Calculated reflectivity curves for an isotropic symmetrical Si400 crystal reflec¬ 
tion at photon energy E = 17479 eV (cr-polaiized) for three curvature radii: 
R = 5.7 m (solid line), R = 2.7 m (dashed line) and R = 1.1 m. The 
inset displays a sketch of the geometty. Results are in agreement with (Erola 
etal, 1990). 


3.2. Laue cylindrically bent crystal in meridional plane 

We calculate here the case described in (Schulze & Chap¬ 
man, 1995) and treated extensively in (Schulze, 1994). It con¬ 
sists of a Silll crystal of thickness T = 700 pm diffracting 
at £ = 33170 eV, bent with meridional Radius R = 325 cm 
(convex to the beam). The crystal directions for ax = 0 are 
n = 1 1 1, Valong = 110, and v_l = 1 1 2. The asymmetrical cut 
is ax = 296.2° (a = —116.2°, x = —206.2°). As discussed by 
Schulze, the anisotropy in the crystal is an important factor to 
be considered. Figure 3 illustrates this by comparing the calcu¬ 
lated diffraction profiles for three cases; i) the isotropic crystal, 
ii) the original crystal configuration (with = 1 1 2), and iii) 
another crystal cut (vj^ = 110). It can be shown here that the 
curves generated by the two methods are in quite good agree¬ 
ment. Fig. 4 shows the variation of the diffraction profile as a 
function of the bending radius. 
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Figure 3 

Diffraction profiles for a meridionally bent {R=325 cm), T — 700 fim thick 
silicon crystal at E — 33170 eV. The calculations show the results for the 
asymmetric {ax — 296.2°) 111 reflection for the isotropic crystal (green lines, 
V — 0.274) and two different crystal cuts: i) vj^ = 112 (red lines) and ii) 
vx = 1 i 0 (blue lines). The solid lines have been calculated using the Penning- 
Polder model and the dashed lines used the multilamellar model. In this case, 
I 77 = 3.Slfirad. The inset displays a sketch of the geometry. 
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Figure 4 

Diffraction profiles computed using the Penning-Polder model for a crystal with 
parameters like in Fig. 3 (vj^ = 1 1 2) as a function of the bending radius. 


The compliance tensor for the chosen crystal cut and asym¬ 
metry angle is shown in Table 1. The variation of the compo¬ 
nents of the compliance tensor versus asymmetry angle that 
affect the diffraction profile are in Fig. 5. 































Table 1 

Compliance tensor values (xlO“^^ rrP-IN) for an asymmetric {olx — 296.2°) 
silicon 111 ciystal. In bold the elements that affect the /3 and c parameters for 
pure meridional bending. The underlined elements contribute to the anticlastic 
bending. 
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Slj 

S2j 

S3j 

■^4; 


^6j 

1 

5.920 

-1.081 

-1.439 

-0.465 

-0.733 

1.489 

2 

-1.081 

6.092 

- 1.611 

1.225 

-1.037 

-0.871 

3 

-1.439 

- 1.611 

6.450 

-0.760 

1.769 

-0.618 

4 

-0.465 

1.225 

-0.760 

14.715 

-1.236 

-2.074 

5 

-0.733 

-1.037 

1.769 

-1.236 

15.404 

-0.930 

6 

1.489 

-0.871 

-0.618 

-2.074 

-0.930 

16.836 


Table 2 

Inputs for a meridionally bent crystal.^? = 10530cm/7o 



Crystal 

Si 111 

Si 311 

Photon energy 

70keV 

Crystal thickness 

0.5 

cm 

Meridional Radius (convex to the beam) 

12759 cm 

10551 cm 

Asymmetry ax 

234° 

263.5° 

Asymmetry x 

1 

0 

-173.5° 


« = 111 

n = 1 1 3 

Crystal cut (for ax — 0) 

^along = 112 

^’along — 3 3 2 


Vx = I 1 0 

vx = I 1 0 


X[deg] 



180 210 240 270 300 330 360 

ax [deg] 


Figure 5 

Variation of the elements of the compliance tensor that affect diffraction for 
meridional bending versus asymmehy angle ax or x- Continuous lines: ele¬ 
ments affecting c and /3. Dotted line: element affecting anticlastic bending. 


3.3. Optimization of a Laue cylindrically bent crystal in merid¬ 
ional plane 

We calculate here the diffraction profile for Laue crystals 
to be used by the high energy X-ray monochromator pro¬ 
posed for the Upgrade ESRF beamline IDS 1 (UPBL02) (ESRF- 
UPBL02, 2012). The monochromator holds two bent Laue 
crystals for monochromatizing the beam in the photon energy 
range 50-150 keV with variable energy resolution in the range 
AEIE « lO”'* — 10“^. controlled by the crystal bending. The 
diffraction plane is horizontal. We simulate a single silicon crys¬ 
tal at energy E = 70 keV using either the reflections 111 
(for low energy and low resolution) or 113 (for high resolu¬ 
tion and high energy applications). The crystals parameters are 
in Table 2. Fig. 6 shows the resulting diffraction profile for the 
111 reflection and Fig. 7 considered crystal configuration. 



E [eV] 


Figure 6 

Diffraction profile for the Sil 11 crystal with parameters in Table 2. 
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Figure 7 

Diffraction profile for the Sil 13 crystal with parameters in Table 2. 


The crystal has been optimized in the following way: i) the 
crystal cut is chosen in such a way that = I 1 0 so the com¬ 
pliance coefficient sgi and are zero for any asymmetry angle. 
Therefore, the crystal is not twisted thus keeping Ic, ™and H 
in the 23 plane. The crystal is cut with an asymmetry angle that 
on one side permits to use both Sil 11 and Sil 13 for low and 














































high resolution applications, respectively; and on the other side 
is optimized for giving good integrated reflectivity. Fig. 8 shows 
the integrated reflectivity and energy bandwidth as a function of 
asymmetry angle and thickness, with indication of the selected 
working points. 




T [cm] 


X[deg] 

-90 -120 -150 -180 -210 -240 -270 
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Figure 9 

Variation of the elements of the compliance tensor that affect diffraction for 
meridional bending versus asymmetry angle ax or x Sit 11. Continuous 
lines: elements affecting c and jS. Dotted line: element affecting anticlastic 
bending. 


3.4. Laue cylindrically bent crystal in sagittal plane 


Figure 8 

Left: Dependency of integrated reflectivity (red) and Ai? (blue) for Sil 11 (solid) 
and Sil 13 (dotted) as a function of the asymmetry angle ax (left). Right: Vari¬ 
ation of integrated reflectivity (red) and Ai? (blue) for Sil 11 at 70 keV (solid) 
and 50 keV (dotted) as a function of crystal thickness (right). 


The same physical crystal should be used for 111 and 113 
reflections, and its compliance tensor is shown in Table 3. The 
variation as a function of the asymmetry (measured for Sil 11) 
of the components contributing to the diffraction is in Fig. 9. 

Table 3 

Compliance tensor values (xlO~*^ rn^/N) for an asymmeh'ic {ax = 234°) 
Sim crystal (or ax = 263.5° Sil 13, corresponding to the same crystal cut). 
In bold the elements that affect the 0 and c parameters for pure meridional 
bending. The underlined elements contribute to the anticlastic bending. 


J 

Slj 

S2j 

^3j 

S^j 

^5; 

S6j 

1 

5.920 

-1.958 

-0.562 

1.071 

0.000 

0.000 

2 

-1.958 

7.010 

- 1.651 

- 1.810 

0.000 

0.000 

3 

-0.562 

- 1.651 

5.613 

0.739 

0.000 

0.000 

4 

1.071 

- 1.810 

0.739 

14.554 

0.000 

0.000 

5 

0.000 

0.000 

0.000 

0.000 

18.914 

2.141 

6 

0.000 

0.000 

0.000 

0.000 

2.141 

13.326 


In this section we analyze a Sil 11 crystal 0.07 cm thick used 
to focus a 50 keV X-ray beam in sagittal direction (Shi, 2011). 
The diffraction plane in the meridional direction is affected by 
the anticlastic curvature originated by the sagittal radius, which 
is much smaller that the typical radii used in meridional focus¬ 
ing. One of the main roles of the crystal is to focus the beam 
in the sagittal direction, and to match the beam divergence in 
the meridional direction (Rowland condition). The optimized 
meridional radius is normally more than ten times larger 
than the sagittal radius required for focusing (e.g., in (Shi 
et al., 2013)). Therefore, in most cases the Poisson’s ratio i/ 
must be minimized (R„, = RJiy, with, Rs = Ri the sagittal 
radius and ly = —si 2 /sii)- Fig. 10 shows the variation of the 
Poisson’s ratio for different crystal cuts. From this graphic it is 
selected the crystal cut (n = I 1 1, Vaiong =211, and Vperp = 
oil) and asymmetry (ax = 125.26° or x = —35.26°). 
The compliance tensor is shown in Table 4. Fig. 11 shows the 
diffraction profile. 
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Figure 10 

Variation of Poisson’s ratio versus asymmetry angle ax (deg) 



Figure 11 

Diffraction profiles for a sagittaly bent Sill I crystal (see text) calculated 
using the Penning-Polder (red) and multilamellai' (blue) models. In this case, 
Ir; = 2.57 firad. 


Table 4 

Compliance tensor values (xl0“*^ nP-IN) for asymmetric (ax = 125.26°) 
Siin crystal. In bold the elements that affect the 0 and c parameters for pure 
meridional bending. The underlined elements contribute to the anticlastic bend¬ 
ing. 


J 

Slj 

Sij 


54; 

S5j 

SQj 

1 

5.920 

- 0.380 

- 2.140 

0.000 

0.000 

0.000 

2 

- 0.380 

5.920 

-2.140 

-0.000 

0.000 

0.000 

3 

- 2.140 

-2.140 

7.680 

0.000 

0.000 

0.000 

4 

0.000 

0.000 

0.000 

12.600 

0.000 

0.000 

5 

0.000 

0.000 

0.000 

0.000 

12.600 

0.000 

6 

0.000 

0.000 

0.000 

0.000 

-0.000 

19.640 


In the last section several diffraction profiles corresponding to 
Bragg and Lane geometries are considered. Numerical calcu¬ 
lations have provided the diffraction profiles. It is possible in 
some cases to obtain some physical parameters like the width 
of the reflections, energy resolution and integrated intensities 
from simple analytical expressions obtained from the approxi¬ 
mated methods applied. Some of these results are discussed in 
this paragraph. 

For Bragg curved crystals the multilamellar is the only 
method that can be applied among the ones described here. The 
shape of the diffraction profile is approximately triangular in the 
case of ’’thick" crystals. The integrated reflectivity 7?^ increases 
with curvature (inverse of radius). The integrated intensity (as 
a function of the p variable) varies from 7?^ = tt for the non¬ 
absorbing thick crystal (77^ = 8/3 for the Darwin solution) to 
the kinematical limit — 7r^|‘FH|/(2A/i) (mosaic crystals). 
For instance, for the Si400 calculated in Fig. 2R^ — 85.47 and 
Rr^ = 20.66,36.44,57.47 for radii R — 5.7,2.7,1.1 m (num¬ 
ber of lamellae 1074,2269,5570), respectively. In case of thin 
crystals, the diffraction profile does not decrease asymptotically 
to zero but decreases abruptly when at a given angle related to 
the crystal thickness. It produces a trapezoidal-shaped profile 
for Bragg crystals, and an almost rectangular shape for Laue 
crystals (as seen in all Laue examples discussed in last section). 
The width of the reflection is (from Eq. 22) Arj = cA = cT A/2. 
The energy bandwidth is calculated using the derivative of the 
Bragg law and Eq. 17; 


T 


AE 

T 


= cot 9bA9 = 


■kPH'Vh? 

2A|7o|sin2 0B^^ 


(31) 


Results of values given by Eq. 31 to the Laue crystals dis¬ 
cussed in the last section are shown in Table 5. 

Eor Laue crystals the PP method gives a diffraction pro¬ 
file width Arj k, fjT in case of low-absorbing crystals curved 
enough to produce a diffraction width much larger than the per¬ 
fect (undistorted) crystal (Arj ^ 2). Using Eq. D.9 it gives a 
bandwidth of: 


A E 2k°7o sin^ 9b 


Table 5 

Comparison of energy bandwidth for different Laue crystal reflections analyzed 
in previous section. 



Fig. 3 

Fig. 6 

Fig. 7 

Fig. 11 

AE FWHM from ML profile 

142.9 

138.0 

30.6 

64.8 

AfiEq. 31 

142.5 

141.6 

31.3 

68.8 

AE FWHM from PP profile 

143.5 

140.6 

30.9 

69.0 

AE Eq. 32 

143.5 

141.3 

31.1 

68.8 


The numerical values for energy bandwidth given by the ana¬ 
lytical formulas agree very well with the calculated ones. In 
fact, both models agree if j3T = cA which implies that 


daz 1 d^{H.u) 
dt 'yok^ dVodVn 


4. Discussion and conclusions 


Therefore, both models are consistent in giving the same AE, 
they give approximated peak reflectivities (checked numerically 































for the treated cases), but will produce different ray trajectories 
(not discussed here). 

We compare now the diffraction profiles produced by the 
same crystal curved along the sagittal direction with the same 
crystal curved with the same radius R =-125 cm along the 
meridional direction, or in both direction (spherical). The 
results are shown in Fig. 12 showing remarkable differences. 
The broadening of the diffraction profile is produced by the cur¬ 
vature in the diffraction plane, therefore the broadening is very 
important for pure meridional or spherical bending. For sagittal 
bending, the broadening is much smaller, because the curvature 
in the meridional plane is due to the larger anticlastic radius 
Rantidastic = —Rmivs, where Rm is the meridional bending and 
Vs = —ii 2 /sii = 0.064 is the Poisson’s ratio in sagittal direc¬ 
tion. In fact, the diffraction profile calculation for sagittally bent 
crystal of R = —125 cm is very well approximated by a merid¬ 
ional bending with Rantidastic = 1947.37 cm (see Fig. 12). 
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Figure 12 

Comparison of the reflectivity curve for the crystal in Fig. 11 (pure sagit¬ 
tal bending, i?., = —125 cm solid line), with the pure meridional bending 
{Rm = —125 cm, dotted line), and spherical bending {Rs = Rm = —125 cm 
dashed line). In red, the diffraction profile produced by a crystal bent in the 
meridional direction v/ith Rantidastic = 1947.37 cm. 


This paper concerns distorted crystals in which the deforma¬ 
tion is created by bending the crystal, usually elastic bending. In 
most cases the crystal curvature is created with the aim of focus¬ 
ing or collimating the X-ray beam, usually in the meridional 
plane, but sometimes in the sagittal direction. We discussed the 
general case of curvature in both directions, which is of inter¬ 
est for spherical bending (equal two-moments), toroidal (dif¬ 
ferent two-moments), or when cylindrical bending is requested 
(one-moment bending). In this last case, however, the elastic 
properties of the materials induce the anticlastic curvature in the 
perpendicular direction. Another important source of deforma¬ 
tion for crystals in synchrotron beamlines is the thermal load. 
Although a full analysis of the crystal reflectivity for heat load 
in crystals is out of the scope of this paper, the same approxi¬ 
mated methods to calculate the diffraction profiles can be used. 
In fact, the original work of (Penning & Polder, 1961a) also 
gives the value of /3 for a crystal deformed by a uniform temper¬ 


ature gradient. Also, our code can be used to estimate if the ther¬ 
mal deformation affects the diffraction properties of the crystal 
by calculating the diffraction profiles for a curvature radius R 
due to the thermal bending of the crystal (Kalus et al, 1973): 


T 

Aj ots 


(34) 


where T is the crystal thickness, ay is the thermal expansion 
coefficient of the crystal and Aj- is the temperature difference 
between the two faces of the crystal that is approximately pro¬ 
portional to the absorber power Aj- = PT/k, with k the diffu¬ 
sion coefficient. This only gives a first estimation of the possible 
alteration of the diffraction profile. Other effects must be con¬ 
sidered, like the thermal expansion of the crystal unit cell, and 
the changes in the rocking curves because of the reflection in a 
non-planar crystal surface (Zhang et al., 2013). 

In conclusion, we summarized the formulation of two 
approximated methods for calculating diffraction profiles: the 
multilamellar applied for both Bragg and Laue bent crystals, 
and the Penning-Polder only applicable to Laue crystals. We 
obtained the general expression of the constant strain param¬ 
eter {f3 for PP and c for ML) including crystal anisotropy for 
any asymmetric crystal cut. The general expressions obtained 
reduce to the formulations of particular cases from literature. 
Last, some approximated expressions for obtaining the angular 
and energy bandwidths of these crystals are given. 

The equations are implemented in the XOP package and are 
validated by studying some examples analyzed in literature. 
The XOP user has access to these calculations using the usual 
XCRYSTAL application. The input files for XCRYSTAL with 
the input parameters of the cases studied in this paper are avail¬ 
able in the examples directory of the XOP distribution. The 
code is also provided in open source under the GPL license at 
https://github.com/srio/CRYSTAL. 

















Appendix A 

Geometry, scattering vectors and angles 

We define an orthonormal reference frame intrinsic to the 
crystal with origin at the crystal center (usually where the cen¬ 
tral ray intercept the crystal) and three vectors {n,Vaiong,v±), 
being n a unity vector normal to the crystal surface pointing 
outside the crystal bulk, and Vaiong and v_l two orthonormal vec¬ 
tors in the plane of the crystal surface, usually chosen to be 
inside and perpendicular to the diffraction plane, respectively. 
A crystal cut is dehned by the asymmetry angle, either a or x- 
Usually a is used in Bragg geometry and x in Laue geometry, 
but if both values are well-dehned they can be used indistinctly 
in both geometries. 

For the formulas used in this text, implemented in the XOP 
code, we have used the following dehnitions and conventions: 

• For simplicity, the three reference vectors previously 
dehned match a simple laboratory reference system cho¬ 
sen n = (0,0,1), Vaiong = (0,1,0) and vj. = (1,0,0). 
The normal of the rehecting crystal surface n is pointing 
outside from the crystal. A generic point in the crystal has 
as coordinates r = (xi,X 2 ,X 3 ). 

• We dehne x as the angle from X 2 axis to H (in mathemat¬ 
ical sense, positive if counterclockwise (ccw)). The angle 
a = X + 90° goes from the crystal surface to the the 
Bragg planes hkl (measured from X2 axis to the rehecting 
surface, thus for symmetric Bragg case a = 180°). The 
asymmetry angle ax dehned in XOP (from crystal sur¬ 
face to crystal planes, positive if clockwise (cw)) holds 
ax + a = 180°. See Fig. 13. 

• The reciprocal lattice vector H can be obtained by rotat¬ 

ing the normal to the surface n an angle e = x- 90 ° = 
—ax around the vector v±. H — (n, 0). The ROT 

operator is implemented via the Rodrigues formula that 
rotates a vector V an angle 9 around an axis a (normal¬ 
ized, 6 is positive in the screw (cw) sense when looking 
in the direction of a) to obtain Vrot'- 

Vro, = ROTs(V, 9) = V cos 0 + 

(a X V)sin9 + a (a ■ V)(l — cos9). (A.l) 

• The (unsigned) Bragg angle verihes A = Idhki sin 0b, 
with A the photon wavelength in vacuum. The direction of 
an incident beam that fulhls Bragg’s law makes an angle 
0B with respect to the Bragg planes, therefore its direction 
is: yo = ROTj^ (-H, 90° - 9b). 

The diffraction prohle for a given hkl rehection is obtained by 
scanning the direction y° of a monochromatic collimated inci¬ 
dent beam (a plane wave) in the vicinity of in the diffraction 
plane. If both vectors are separated by an angle A0 = |0| — \0b\ 
we can set = ROT^^(Vb, —A9). 

The vector expressions can be expressed in angles, but care 
must be taken with the definition and sign of angles. 

The incident angle (measured ccw from X2 to V^) for the ray 
fulfilling the Bragg law (in both Bragg and Laue geometries) is 
00 = 180° + a — 0B, and the reflected angle (measured ccw 


from X 2 to V^) is 9h = 180° + a — 0 b. In Laue geometry, for 
a given crystal cut, we may send the beam to match the Bragg 
angle below or onto the Bragg planes (which strictly speaking 
correspond to using ztH), and reversing time. From these four 
cases only two are independent. 



Figure 13 

Schematic view of diffraction using asymmetric Silll crystals at £ = 

4 keV, ds = 29.62°. Left: Laue geometry: ax = 125.26°, x = —35.26°, a = 
54.74°. Right: Bragg geometry: ax = 10°, x = 80°, a = 170°. The plotted 
incident direction comesponds to the Bragg position Vj), so here =12- 


The directions (as as function of a, ax = 180° — a, and 
X = a — 90°), for the directions fulfilling the Bragg law (Eq. 1) 
are: 











= (O,cos6»o,sin6»o) = 

(0, — cos(0B — a:), sin(0B — a:)) = 
(O,cos(0B + ax), - sin(dB + ax)) = 
(0, - sin(6»B - x), - cos(0B - x)) 

Vb = (0,cos6*H,sin6»H) = 

(0, - cos(0B + a), - sin(0B + a)) = 
(O,cos(0B - ax),sin(6»B - ax)) = 
(0, sin(6»B + x), - cos(6»b + x)) 

= (0, sin a, — cos a) = 

(0, sin ax, cos ax) = 
(0,cosx,sinx) 


with M, the displacements (elongations) and x, the coordinates 
along the three spatial dimensions i = 1,2,3. 

Then, Eq. B.2 can be expanded as: 


(A.2) 


du2 
3x3 
3 m 1 
3x3 
3 m 1 
3X2 


3mi 

3xi 

3m2 

3x2 

3m3 

3x3 

3m3 

3x2 

3m3 

3xi 

3m2 

3xi 


■Silo’ll + ■S120’22 + ■S130’33 + .Sl40’23 + ■Sl50’53 + .S160’52 

■S2lO’ii + .S220’22 + *23033 + .S24O23 + ■525013 + ■S26O12 

■S3IO11 + ■S32O22 + ^533033 + S 34 (J 23 ■S35O13 + ■S36012 

■S4ICTjj + ^542022 + ■S43O33 + S 44 (J 23 ■S45O13 + ■S 46 oj 2 (B. 5 ) 

■S5IO11 + ■S52O22 + ■S53O33 + S 54 (J 23 ■S55053 + ■S 56 Crj 2 

■S6IO11 + ■S62O22 + ■S63O33 + ■S64O23 + ■S65O13 + ■566012 


Appendix B 

Elastically bent crystals 


B.l. Bending an anisotropic plate 

The generalized Hooke’s law makes a “linear” relation 
between the the strain tensor (5'y, i,j = 1,2,3, characteriz¬ 
ing the deformation, adimensional), and the stress tensor 
k,l = 1,2,3, generalized forces with dimension of force times 
L-2) (Eq. 1.3.2 in (Hearmon, 1961)): 

~ ^'ijkl^'kl (B-1) 

where is the elastic compliance tensor. The stiffness tensor 
is the inverse of the compliance tensor, but here we will center 
the discussion on the compliance tensor only. Each index can 
take three values (1,2,3, corresponding to the three directions 
in 3D space) and the repeated indices are summed. Erom sym¬ 
metry considerations, not all 81 elements of the compliance are 
independent, but they are reduced to 36. Moreover, thermody¬ 
namical considerations reduce the number of independent com¬ 
pliance elements to a maximum of 21 in the most general case. 
These facts make possible to reduce Eq. (B.l) to a simpler form 
(Eq. 1.3.6 in (Hearmon, 1961)): 

Sq = SgrO-r {q,r = 1,...,6), (B.2) 

with now Sqr a 36 elements symmetric matrix. The new strain 
6-dim vector as a function of the old strain tensor is (Eq. 1.2.8 
in (Hearmon, 1961)): 


‘^1 ~ ‘^lli ‘^2 — ‘^221 ‘^3 — ■^ 33 ) 

54 = 25^3, 55 = 25;3 , 56 = 25'i2, (B.3) 


(the same convention applies for indices from s' to s, e.g., 
*36 = ■ 53312)5 (Eq. 1.2.4 in (Hearmon, 1961)) 


‘’j 2\dxj dxi) 


(B.4) 


In general, the stress tensor is proportional to the generalized 
torques M (per unit of length) and forces A. Let us suppose a 
crystal as a rectangular plate of thickness T with axes 1 and 2 
parallel to the the edges, and axis 3 perpendicular to the surface. 
Considering only two torques Mi and M 2 (dimensions torque 
per length), applied in-plane (Eig. 14), we have: 


, M 1 X 3 

, M 2 X 3 

^22 — — j — 

(T 33 = 0 

^13 = ^23 = ‘^12 = 0 


(B.6) 




X 2 




Figure 14 

Double moment bending of a rectangular plate of thickness T 



Replacing all into Eq. (B.5) we obtain; 

dut 1 , , 

— = - (S11M1X3 + S12M2X3) 

0 X 1 / 

3m2 1 

— = - (S21M1X3 + S22M2X3) 

0X2 1 

^ = 7 {s2\M\X2 + S32M2X3) (B.7) 

0X3 / 

du2 3m3 1 

7 - 1 - 7— = 7 (* 41 ^ 1 X 3 + S42M2X2) 

0X3 0X2 I 

dui dm 1 ^ ^ 

7- 1 - 3— = 7 (* 51 ^ 1 X 3 + S52M2X2 ) 

0X3 OX 1 1 

dui du2 1 / \ 

7 -b 7 — = 7 (■S' 61 -^ 14 i ;3 + S 62 M 2 X 2 ) 

0x2 <jx 1 1 

Integrating these equation we obtain (Eq. 4 in (Chukhovskii 
etai, 1994)): 


Ml = +Si2M2)xiX3 

+ (s5iMi + S52M2)x^/2 + {saiMi + 5624^2 )x2X3/2] 

M2 = y[(s2lMi +S22M2 )x2X3 

+ (^41471 + S42^2)x^I'^ + (sei-^l + ■S'624^2 )''CiX3/2] (B. 8 ) 

M3 = ^[—(sil-^l + ■S'l2-^2)4ti — (■S'21-^1 + ■S'22-^2)-t;2 

-(s6i47i + 562^2 )xiX2 + (s3lA^l + S22M2)xl] 

Replacing X 3 = — 7’/2 we obtain the equation of the plate 
surface (Eq. 5a in (Chukhovskii et al., 1994))): 


B.2. The compliance tensor for crystals cut along given direc¬ 
tions 

The number of independent components in the compliance 
tensor depends on the crystal symmetry. Eor a crystal cut along 
the crystallographic axes, triclinic crystals have all 21 compo¬ 
nents independent, monoclinic crystals have 13, orthorhombic 
crystal 9, tetragonal crystals 6 or 7 and cubic crystals 3 (see 
table 5 in (Hearmon, 1961)). An isotropic material has only two 
independent components. Eor a crystal cut along given direc¬ 
tions, fourth-rank compliance tensor must be transformed fol¬ 
lowing (Eq. 1.5.2 in (Hearmon, 1961)): 

^ijkl ~ ‘^imOjnOkoOlp^mnop (B.12) 

where the elements of the a tensor are the direction cosines of 
the new axes. Software codes are available (Honkimaki, 2014) 
(Schulze & Chapman, 1995) to transform a generic compli¬ 
ance tensor expressed in the reference frame coincident with the 
crystallographic axes, to a new one where the crystal has been 
cut following known planes expressed by their Miller indices. 

Eor a cubic crystal cut along its crystallographic axes, the 
compliance tensor has only three different values. The non¬ 
zero elements are: S12 = ^21 = sb = ssi = ■S23 = ■S32, 
■sii = *22 = *33 and S 44 = S 55 = S 66 - Because of cubic sym¬ 
metry, this tensor is independent on the choice of the axes in 
the crystal (in plane or normal). The components of the compli¬ 
ance tensor for the most usual crystals are expressed in Table 1. 
The components of the compliance tensor for a cubic crystal cut 
along a given direction can be calculated using the generic Eq. 
(B.12), resulting in analytical expressions given by (Wortman 
& Evans, 1965). A convenient form easy to implement in com¬ 
puter languages is given in (Zhang, 2010) (Zhang et al, 2014) 


x^ x? 

M3 = — (silMi +512^2)^ — {s2\Mi + S22M2) 2^ 

— (sfiiMi +562^2)-^ + (ssiMi + S32M2)— (B.9) 

ZI ol 

Erom here, one can obtain the prohles along the i = 1,2 
directions, which are approximately circular with the form M 3 = 
—x?/(2R,). The radii are: (Eq. 5b in (Chukhovskii et al., 1994))): 


1 Ml M 2 

— — ‘^11- “b - 

Ri I I 

1 Ml M 2 

Sn—+S22 — 


(B.IO) 


Erom Eq. (B.IO) is possible to calculate the applied torques 
as a function of the curvatures: 


Ml 

1 1 

f «12 

* 22 ' 

/ 

S 12 S 21 — S 11 S 22 

\R2 


M 2 

1 

( S 2 I 

ill 

/ 

S 12 S 21 — S 11 S 22 

Ul 

R 2 


(B.ll) 


The Eq. B .8 give the displacements m, knowing the elastic 
compliance tensor and the torques applied (or more interest¬ 
ingly, via the bending radii using Eq. B.IO). 


In the particular case of one single torque (M 2 = 0) Eqs. 
(B.IO) give the main bending radius R\ = II{s\\Mi) and a cur¬ 
vature radius in the perpendicular direction R 2 ~ Ri(siiAi 2 ), 
the anticlastic curvature. The ratio of the transverse component 
of the compliance tensor (in this case along the direction 12 ) 
over the longitudinal component is the Poisson’s ratio v \2 = 
—S 12 A 11 . Eor example, a silicon crystal cut along the crystal¬ 
lographic axes has (see table 1) v \2 = — (—2.14)/7.68 = 0.28 
thus the anticlastic radius is 3.6 times larger than the main bend¬ 
ing radius with the curvature in opposite direction. Note that the 
Poisson’s ratio and the anticlastic radius in a crystal cut along 
directions different from the crystallographic axes are different. 

Table 1 

Compliance tensor i elements (x 10“*^ m^lN) for most used cubic perfect crys¬ 
tals (Si, Ge: (Wortman & Evans, 1965), Diamond: (Berman, 1965); see review 
in (Hedayat et a/., 2012)). For completeness, the relationships between the com¬ 
pliance i and the stiffness c tensors for cubic crystals are included. 



Si 

Ge 

Diamond 

s and c relationships 

^11 

7.68 

9.64 

1.04 

(cii -bci 2 )/[(cii -Ci 2 )(cn + 2 ci 2 )] 

‘?12 

-2.14 

-2.60 

- 0.211 

-Cl 2 /[(cil -ci 2 )(cii + 2 ci 2 )] 

S44 

12.6 

14.9 

1.93 

I/C 44 

Cl 2 




-il 2 /[(-Sll -■« 12 )(ill + 2 ^ 12 )] 

cn 




(ill + il 2 )/[(ill - il 2 )(ill + 2 ^ 12 )] 

C 44 




l/i44 



Table 2 


Equivalence for the compliance tensor indices for different reference frames 
used in literature:Shi (Shi, 2011), Honkimaki (Honkim^i, 2014), Schulze 
(Schulze, 1994) (Schulze & Chapman, 1995), and Zhang (Zhang et al, 2014). 


This Paper 

Shi 

Honkimaki 

Schulze 

Zhang 

1 

3 

1 

2 

2 

2 

-2 

-3 

1 

3 

3 

1 

2 

-3 

1 

4 

6 

4 

5 

5 

5 

5 

6 

4 

6 

6 

4 

5 

6 

4 


Appendix C 

The parameter c in the multilamellar model. 

In the multilamellar model, for diffraction occurring in the 
X 2 ,X 3 plane, the bent crystal is approximated by a stack of 
lamellae which gradually change direction and lattice spacing. 
The reduced curvature c can be expressed as: 


dn dri daz dt bA daz ,, 

c = — =-= —;= -(C.l) 

dA daz dt dA 

where the definitions of the 77 and A in (Zachariasen, 1967) have 
been used, and f = X 3 . 

Taupin (Taupin, 1964c) discussed the variation of az as a 
function of the crystal deformation and form of the incident 
wave, which can be chosen in such a way that az is only depen¬ 
dent on the thickness direction t = X 3 . In this case, he found a 
tensorial expression ((Taupin, 1964c) Eq. II.1.7): 


da^ 

dt 


2 Ml M 2 ,. , 

-[(*21— + *22 —Mi + 

70 f 1 


Ml ■^2 n, , . Afr M 2 , , 

(sai — -f S32—j-)A2 + (■S'41 — -1" Si\2^-)A2,\ 


(C.4) 


with: 


Al = -vMiy? - y^) + y^y^ivi - ^2°)+ yMiy? - y!) 

(70 - 7h)(1 + 7o7h) 

A 2 = y^y^{y^ - y 3 '’) = 7o7/r(7H - 70 ) 

A 3 = y^y^{y? - y 2 “) = 

1o1h{\I 1 - 7i - \/l - 7 o) (C-5) 

For computing A 3 it is preferred to use the form without the 
square toot, to avoid incertitude due to the double sign. For par¬ 
ticular case widely treated in literature of meridionally bend¬ 
ing Ml = 0, M 2 II = l/(f? 2 'S 22 ) and isotropic (Poisson’s ratio 
V = —S 23 /S 22 , and S 42 = 0 ) we have: 


da. 


-[Al - iyA2] 


dt 7o 7?2 
which gives a c parameter like (Caciuffo et al, 1990): 


(C. 6 ) 


A(b- l)|7tf| 


[1 +b{l +v)ll] , 


For the particular case of sagittal bending, M 2 
Mill = l/(/?iSii) the Eq. C.4 reduces to: 


(C.7) 
0 and 


daz 

dt 


oxjdxk 

K—1 J—l (—1 


(C.2) 


daz 

dt 


2 .S21 , , ■531 , , *41 , X 

—Al -f —A2 -f —A3) 

7o«i ill ill ill 


therefore: 


(C. 8 ) 


where y°’^ = y 3 ®'^) are unitary vectors along the 

incident and diffraction directions, as defined in the text. 

Performing the summation in the diffraction plane 23 (thus 
y^ = y° = 0 ) one gets: 


daz 

dt 




70 


V2X« - 

y3X(y2" 

v^y^{y^ 

y^yiiyi 


V«)^+ 


■y| 


y| 


■V2^) 


3x23x3 

. 3^m2 

3 X 3 

3^ 
3 x 3 
03 3^«2 


-f 


-f 


3x23x3 ' 


(C.3) 


Inserting the derivatives calculated from the equations for the 
displacements B. 8 , one obtains: 


I'A^I ^ ^"21 


S 3 I , , ■*41 , N 

-A2 -f -A3) 

ill ill 


(C.9) 


Appendix D 

The parameter j3 in the PP model 


The strain gradient (3 of an asymmetric Faue crystal with 
respect to the crystal surface is defined by 

where G is defined as: 


Q ^ d^jH.u) 

~ dyodyn 


(D.2) 


For calculating G we restrict the calculation to the diffrac¬ 
tion plane 23, thus replacing H.u = (n^M 2 + n 2 U 2 ,)/dhki, To = 
( 0 , y)*, y 3 ‘’) and Vh = ( 0 , y^^y^ ) we obtain: 



2 sin X tan 6b 1 


G={%.Vr)(VH.Vr){H.u) = 


O’l/HH 


dhkl 




M2 


Ml 


*23 H-pSl3 I 


M2 


(,^2 ^3 ”2 + ^3 ^2 ”2 ~ ^2 ^2 ”3 j 1 —^22 + 

(M2 Ml \ 

[(jl I —S 22 + —^13 I + (j2 

M 2 Ml 


/ 

Ml 


*24 ■ 


1 


Af2 

—^5'24 H- 


G 3 


-*22 ■ 


-*12 


Ml 


-*14 


*12 = 


(D.3) 


0 H 

Considering that 70 .// =V^’ , and replacing the vector com¬ 
ponents by their angular expressions (Eq. A) we get: 


Gi = = 7o7h sin X 

G 2 = = 7 o7h cos x 

^ _ t/Ot/H H I t/Ot/H H \ 7 O\ 7 H H _ 

G3 — V2 ^3 ^2 + ^3 ^2 ^2 ~ ^2 ^2 ^3 “ 


-(1 +707//) sinX = - ( 1 + 


cos 26b + cos 2 x 


(D.4) 


sinx 


For bending the crystal only in meridional direction with a 
curvature radius Rm = R 2 we have M\ — Q and M 2 II = 
(^ 2 * 22 )”^ (from Eq. (B.IO)). Inserting these values into Eq. 
(D.3) we obtain (note also that 2707 // = cos 26b + cos 2x): 


G = -—— sin V 

dhklR2 


1 + 


cos 26b + cos 2 x 


1 - 


*23 + cot X*24 


*22 


which inserted in (D.9) gives: 


(D.5) 


1 + 


^ 2 sinx tan 1 

cos 26b + cos 2x ( *23 + cot x *24 


1 - 


*22 


(D.6) 


The (3 value in Eq. (D.6) was first obtained by (Schulze & 
Chapman, 1995). The equations are identical considering the 
different choice of the axis. Noticeably, there are typos in the 
indices of the compliance tensor in Eq. 5 of (Schulze & Chap¬ 
man, 1995) (as well as in Eq. (37) of (Schulze, 1994)). If using 
their reference system as defined in Fig. 1 of (Schulze & Chap¬ 
man, 1995) (or Fig. 7 in (Schulze, 1994)) one should obtain 
an equation like Eq. (D.6) but replacing S 22 —* 11 , *23 — 
* 13 , *24 —*15 (see Table 2). The equation as written in these 
two references correspond to a different reference frame, like 
the one shown in Fig. 10 of (Schulze, 1994). 

The case of isotropic materials can be obtained setting S 24 = 
0 in Eq. (D.6) and considering the definition of Poisson’s ratio, 
as minus the ratio of the compliance term element along the 
direction transversal to the curvature and the one along the 
curved direction (i.e., y = —* 23 /* 22 - We obtain (Sanchez del 
Rio et al, 1997): 


P = - 


1 + 


cos 26b + cos 2x 


(1+^) 


^ (D.7) 

The case of single bending in sagittal direction was studied 
in (Shi, 2011). Setting Rs = R\, M 2 = 0 and M\II = (Risn)“' 
(from Eq. (B.IO)) one obtains: 




*12 


-1- 7o7// — 


*11 


*12 *13+COtX*14 


*11 


*11 


which inserted in (D.9) gives: 

2 sin 6b 


/? = - 


1 


-707//- 


Sl3 +COtX*14 


Sin X X 


sn 


+ (1 +707//) — 
*11 


which corresponds to Eq. 18 in (Shi, 2011). 


(D.8) 


(D.9) 
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